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' Abstract. The discovery of the microquasar LS 5039 well within the 95% conf- 

icence contour of the Unidentified EGRET Source (UES) 3EG J1824-1514 was a 
major step towards the possible association between microquasars (MQs) and UESs. 
The recent discovery of precessing relativistic radio jets in LS I +61 303, a source 
' ^ . associated for long time with 2CG 135+01 and with the UES 3EG J0241+6103, has 

' given further support to this idea. Finally, the very recently proposed association be- 

^ \ tween the microquasar candidate AX J1639. 0-4642 and the UES 3EG J1639-4702 

■ points towards a population of High Mass X-ray Binary (HMXB) /Neutron Star 

OO ' (NS) microquasars as counterparts of low-latitude unidentified EGRET sources. 

o : 
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1. Introduction 



o 

h . 

bo . The third EGRET catalog (Hartman et al. 1999) contains 271 point 

I sources detected at energies above 100 MeV. The majority of these 

sources, ~168 or ~62%, still remain unidentified. Among them, there 
are 72 sources located at low galactic latitudes, having |6|<10°, which 
• represents around 45% of the UES population. Therefore, several of 

these objects are presumably of galactic nature. Similar properties be- 
tween some of these UESs, indicate that there are at least three different 
groups of galactic populations (Romero et al. 2004, Grenier 2004). The 
group of young stellar objects and star- forming regions (Romero 2001), 
those sources forming a halo around the galactic center and a group of 
sources correlated with the Gould Belt (Grenier 2000). 

Based both on multiwavelength observations and theory, microquasars 
(see Mirabel & Rodriguez 1999 for a review) with massive companions 
have been proposed as possible counterparts of the first group of galac- 
tic UESs by several authors (Paredes et al. 2000, Kaufman Bernado 
et al. 2002, Romero et al. 2004, Bosch-Ramon et al. 2004). In Sects. 2 
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and 3 of this paper we will briefly review the properties of the two well- 
known microquasars LS 5039 and LS I +61 303, typically associated 
with the first group of UESs, while in Sect. 4 we will present the possible 
association between the microquasar candidate AX J1639.0— 4642 and 
the UES 3EG J1639-4702. Finally, in Sect. 5 we will compare the 
available data of these 3 sources from radio to gamma-rays, and we 
will discuss on similarities pointing towards a population of HMXB 
with NS microquasars as counterparts of low-latitude UESs. 



2. LS 5039 

The high mass X-ray binary system LS 5039 (Parcdcs ct al. 2000) is 
one of the ~15 confirmed galactic microquasars (Ribo 2004). LS 5039 
is a bright V^ll.2 star with an ON6.5V((f)) spectral type (McSwain 
et al. 2004). The binary system has a short orbital period of P = 
4.4267 ±0.0005 d, a high eccentricity of e = 0.48 ±0.06, and a low mass 
function f{m) = 0.0017 ± 0.0005 Mq, suggesting the presence of a NS 
as the compact object in the system (McSwain et al. 2004). 

Observations conducted with the EVN and MERLIN (see Fig. 1) 
confirmed the persistent nature of this MQ, and revealed the presence 
of an asymmetric two-sided jet reaching up to 1000 AU on the longest 
jet arm (Paredes et al. 2002). These observations also suggest a bending 
of the jets with increasing distance from the core and/or precession. 

The possibility that LS 5039 is a 7-ray emitter was suggested by 
Paredes et al. (2000), who proposed the association of the system 
with the UES 3EG J1824-1514 (Hartman et al. 1999). We show in 
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Figure 2. Location map of 3EG J1824— 1514. The contours represent, from inside 
to outside, the 50%, 68%, 95%, and 99% statistical probability that a 7-ray source 
lies within the given contour. The only source with X-ray and radio emission (filled 
circle inside an open big circle, I = 16.88° and b = —1.29°), and well inside the 95% 
contour, is LS 5039 (Ribo 2002). 



Fig. 2 the location map of the 7-ray source together with the NVSS 
and bright /faint ROSAT sources. The only simultaneous X-ray/radio 
source within the statistical contours of 3EG J1824— 1514 is the mi- 
croquasar LS 5039. We note that this binary system is present in the 
BATSE Earth occultation catalog of low-energy gamma-ray sources 
(Harmon et al. 2004), with a positive detection of a few mCrab up 
to ~100 keV. The source is not present in cumulative observations 
conducted with the INTEGRAL satellite (Bird et al. 2004), although 
it is expected to be detected when adding a few more months of data. 
We also point out that there is an unidentified COMPTEL source with 
a position compatible with LS 5039 (Collmar 2004). 

Astrometric studies carried out by Ribo et al. (2002) , show that it is 
a runaway system with a systemic velocity of ~150 km s~^ that moves 
away from the Galactic plane with a velocity of --^100 km s^^. This 
result, combined with the possible lifetime of the donor star, indicates 
that it could reach a not-so-low galactic latitude of 6 = —12° still 
behaving as a microquasar. 

Bosch-Ramon &; Paredes (2004a) have recently developed a detailed 
numerical model to test whether this system can actually produce the 
emission detected by EGRET through inverse Compton (IC) scatter- 
ing. Their numerical approach considers a population of relativistic 
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Figure 3. Location map of 3EG J0241+6103. The position of LS I +61 303 is marked 
with an asterisk. Adapted from Hartman et al. (1999). 



electrons entrained in a cylindrical inhomogeneous jet, which interact 

with both the radiation and the magnetic fields. The computed spec- 
trum is able to reproduce the observed spectral characteristics at very 
high (GeV) energies. 



3. LS I +61 303 

The Be/X-ray binary system LS I +61 303 is a well-studied object since 

it presents radio and X-ray variability linked to its ^^26. 5 d orbital 
period (Gregory 2002; Paredes et al. 1997). The donor star in this 
system is a rapidly rotating BOV star with variable mass loss (Hutchings 
&; Crampton 1981). Some properties of this system can be explained 
assuming that the unseen companion is a non-accreting young pulsar 
with a relativistic wind strongly interacting with the wind of the Be 
star (Maraschi &i Treves 1981). On the contrary, other properties of 
LS I +61 303 fit better a model where the companion is accreting even 
with two episodes of super-critical accretion along the orbit (Marti &; 
Paredes 1995). 

This X-ray binary system has been associated for long time with 
the gamma-ray source 2CG 135+01/3EG J0241+6103 (see Fig. 3), 
which displays variability on timescales of days (Tavani et al. 1996, 
1998; Wallace et al. 2000). During the last years, Massi et al. (2001, 
2004) have revealed its MQ nature through the discovery of a radio 
jet (see Fig. 4) extending 200 AU at both sides of a central core. 
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Figure 4- a) MERLIN self-calibrated image of LS I +61 303 at 5 GHz and using 
natural weights, obtained on 2001 April 22. The synthesized beam is plotted in the 
lower left corner. The S-shaped morphology strongly recalls the precessing jet of 
SS 433, whose simulated radio emission (Fig. 6b in Hjellming & Johnston 1988, 
rotated here for comparison purposes) is given in the small box. b) Same as before 
but for the April 23 run and using uniform weights for better displaying. From Massi 
et al. (2004). 

that appears to experience a fast precession, which could explain the 
short-term gamma-ray variability of 3EG J0241+6103 (as proposed 
by Kaufman Bernado et al. 2002) and the puzzling VLBI structures 
found in previous observations. This result points to the occurrence of 
accretion/ejection processes in this system, ruling out, in principle, the 
non-accreting young pulsar scenario. 

Massi (2004) has recently studied the data acquired within the 
pointed EGRET observations of 3EG J0241+6103 and claimed the 
detection of a periodicity of P = 27.4ib7.2 d, consistent with the orbital 
period of the binary system. If this is confirmed, the identification of 
LS I +61 303 as the counterpart of 3EG J0241-I-6103 would be unam- 
biguous. In any case, an important point is that all the available data 
are compatible with an increase of 7-ray emission around periastron, 
that can be tested with Cherenkov telescopes and future satellites. 

This binary system is also present in the BATSE Earth occultation 
catalog of low-energy gamma-ray sources (Harmon et al. 2004), with 
a positive detection of a few mCrab up to ~100 keV, although the 
detection is not as significant as in the case of LS 5039. LS I +61 303 is 
not present in cumulative observations conducted with the INTEGRAL 
satellite (Bird et al. 2004), although it is expected to be detected when 
adding a few more months of data, like in the case of LS 5039. We note 
that there is a COMPTEL source containing LS I +61 303 and the 
quasar QSO 0241+622 (van Dijk et al. 1996). 
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Figure 5. Contour image of the MGPS data obtained with MOST at 843 MHz 
on 1992 April 16. The radio source MOST J1639. 0-4642 is weU within the 90% 
uncertainty error circle of the X-ray source AX J1639. 0—4642. The eUipse in the top 
right corner is the convolving beam. Prom Combi et al. (2004). 

A numerical model to explain the EGRET emission of LS I +61 303 
has also been developed by Bosch-Ramon & Paredes (2004b). 



4. AX J1639.0-4642 

Aimed at discovering new MQs, Combi et al. (2004) have recently 
carried out a multiwavelength study of the unidentified X-ray source 
AX J1639.0-4642. This object was discovered by the Advanced Satel- 
lite for Cosmology and Astrophysics (ASCA) observatory at the 0.7- 
10 keV energy range, and presented as a possible HMXB (Sugizaki et al. 
2001). Its measured flux was i*x(o.7-io keV)= 19-2 x 10~^^ erg cm~^ s~^, 
and it showed variable X-ray emission, with a confidence > 99 %. Its 
spectrum was fitted with a power law with a very hard photon index 
r = -0.01 

-aeo ^"^^ ^ poorly constrained hydrogen column density of 
Nil = I2.82tlil X 10^^ cm~2_ Combi et al. (2004) re-analyzed these 
data and found evidences for variability on timescales of hours. 

In searching for radio sources in the field of AX J1639. 0—4642, 
Combi et al. (2004) found that the Molonglo Galactic Plane Survey 
(MGPS) at 843 MHz (Green et al. 1999) revealed a point-like radio 
source (see Fig. 5), dubbed MOST J1639.0-4642, weh within the error 
box of the X-ray source, with a fiux density of 136 =b 18 mJy. At near 
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Figure 6. 2MASS /fs-band image of the environment of the radio source 
MOST J1639. 0—4642, whose Sct position error circle is shown, together with the 
2(7 position error eUipse of the far infrared source IRAS 16353—4636. Several NIR 
sources from the 2MASS catalog are possible counterparts of both sources. Prom 
Combi et al. (2004). 

infrared (NIR) wavelengths Combi et al. (2004) inspected the 2 Micron 
All Sky Survey (2MASS, Cutri et al. 2003), and found 10 sources in 
the 3(T error circle in position of MOST J1639. 0—4642, some of them 
visible in the Kg-hanA image shown in Fig. 6. At the far infrared 
part of the spectrum, from 12 to 100 microns, they found that the 
source IRAS 16353—4636 lies inside the error box of the X-ray source. 
This source overlaps the southern part of the 3(J position error circle 
of MOST J1639. 0—4642, and its location uncertainty ellipse contains 
several 2MASS sources, as can be seen in Fig. 6. The X-ray source 
AX J1639. 0—4642 has been recently rc-discovcrcd at higher energies 
with the IBIS telescope onboard the INTEGRAL satellite, dubbed 
IGR J16393-4643 (Malizia et al. 2004; Bird et al. 2004). This source 
shows an average flux of -Fx(20-ioo keV) — 5 x 10^^^ erg cm~^ s~^, and 
presents a factor of 2-3 flux variability on timescales of months. 

Although there is no spectroscopic/photometric information of a 
NIR/optical counterpart to derive a distance to the source, assuming 
that it is located in the Scutum-Crux or in the Norma spiral arms, a 
range of distances between 3 and 13 kpc is obtained. 

Combi et al. (2004) pointed out that AX J1639.0-4642 lies inside 
the 95% location contour of the UES 3EG J1639-4702 (Hartman et al. 
1999) as can be seen in Fig. 7. Its 7-ray flux is (53.2 ± 8.7) x 10~* pho- 
ton cm~^ s~^, presents a steep 7-ray spectral index of F = 2.5 ± 
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Figure 7. Location map of 3EG J1639— 4702. The position of tiio X-ray source 
AX J1639. 0—4642, well within the 95% probability contour, is maxked with an aster- 
isk. Other sources appearing outside the 99% probability contour are also indicated. 
Prom Combi et al. (2004). 

0.18 and has a variability index of I = 1.95. Although Torres et al. 
(2001a) found three radio pulsars inside the 95% confidence contour 
of the 7-ray source, its possible variability and steep photon index 
do not seem to agree, in principle, with a pulsar origin. Similarly, 
these properties would rule out an association with the three SNRs 
found within the 95% confidence contour (Torres et al. 2003). More- 
over, no identified blazar has been found within the 7-ray contours. 
Therefore, Combi et al. (2004) suggested that the microquasar candi- 
date AX J1639.0-4642/MOST J1639.0-4642 (=IGR J16393-4643) is 
the counterpart of 3EG J1639-4702. Observations with ATCA are in 
progress to unveil the nature of this source. 

5. HMXB/NS microquasars as counterparts of low-latitude 
unidentified EGRET sources 

As discussed above, the possibility of MQs being 7-ray emitters was 
suggested by Paredes et al. (2000), who proposed the association be- 
tween the HMXB LS 5039 and the UES 3EG J1824-1514. In their 
scenario (Paredes et al. 2000, 2002) the 7-rays are produced by IC 
upscattering of stellar ultraviolet (UV) photons by the non-thermal 
relativistic electron population that later on will produce the detected 
radio emission. Recently, more detailed models considering precession 
(Kaufman Bernado et al. 2002), hadronic jets in windy microquasars 
(Romero et al. 2003) and all possible photon fields (Bosch-Ramon 
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Table I. Properties of the three 7-ray sources discussed 
in the text. The luminosity interval for AX J1639.0-4642 
reflects the range of assumed possible distances. 



7-ray source 






-'^tOIOO MeV)'' 


d 








(erg s~'^) 


(kpc) 


3BG J1824-1614 


2.19±0.18 


3.00 


3.6x10^*^ 


2.9"* 


3EG J0241+6103 


2.21±0.07 


1.31 


3.1x10^*5 


2.0 = 


3EG ,71639-4702 


2.50±0.18 


1.95 


(3-50) X 10^^ 


3-13 ? 



^ Hartman ct al. (1999); ^ Torres et al. (2001b); computed using the 
photon fluxes and indices from Hartman et al. (1999); ^ Ribo et al. (2002); 
" Frail & Hjellming (1991). 



Table II. Proposed X-ray/optical/radio counterparts of the three 7-ray sources of 
Table I. In the cases of the microquasars LS 5039 and LS I +61 303 the luminosity 
intervals correspond to intrinsic variability of the sources at the corresponding 
distances, while in the case of AX J1639. 0—4642 they reflect the range of assumed 
possible distances. 



X-ray source 


-^X(0.7-10 keV) 
(erg s^^) 


-^radio(n-l-100 GHz) 
(erg s^^) 


Spectral 

type 


(days) 


d 

(kpc) 


LS 5039 


(0.5-5) X 10=** " 


~1.0X10^^ 


ON6.5V((f)) " 


4.4"^ 


2.9'' 


LS I +61 303 


(1-6) XIO-'** ■= 


(1-17)X10-''1 f 


BOVe" 


26.5 


2.0 * 


AX J1639. 0-4642 


(2-40) X 10^* 


(0.8-16) X 10^^ ^ 


? 


? 


3-13 ? 



" Reig et al. (2003); Marti' et al. (1998), Ribo et al. (1999); ■= McSwain et al. (2004); Rib6 et al. (2002); 
= Paredes et al. (1997); ^ computed from the values quoted in Strickman et al. (1998); ^ Hutchings &: 
Crampton (1981); Gregory (2002); » Frail & Hjellming (1991); ^ computed assuming a = -1. 



et al. 2004) have been proposed to explain the high-energy gamma-ray 
emission from HMXB microquasars. 

On the other hand, as already stated, the X-ray binary system 
LS I +61 303 has been associated with the UES 3EG J0241-F6103, 
and Massi et al. (2001, 2004) have revealed its MQ nature. If the MQ 
nature of AX J1639.0-4642 is confirmed, it could be the third MQ 
source related to a UES. 

We quote the basic properties of these three 7-ray sources in Table I, 
and the properties of the proposed X-ray counterparts in Table II. 
There are 3 observational facts that should be noted. The first one 
is that LS 5039 and LS I -^61 303, and probably AX J1639.0-4642, 
have massive optical companions, which provide an intense stellar UV 
photon field. On the other hand, the compact object appears to be com- 
patible with a neutron star in the cases of LS 5039 (McSwain et al. 2004) 
and LS I -|-61 303 (Hutchings & Crampton 1981; Casares et al, 2004; 
but see Massi 2004) . Finally, it is interesting to point out that the lumi- 
nosities obtained in each spectral domain are very similar in all three 
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sources, specially for the shorter distances to AX J1639. 0—4642, giving 
support to the idea that all of them have similar emission processes. 

Moreover, it should be noted that LS 5039 and LS I +61 303 are the 
only microquasars having both a high-mass donor and (possibly) a NS 
as the compact object, and that there are no other microquasars (con- 
taining low-mass donors and/or black holes) located within the prob- 
ability contours of unidentified EGRET sources. Therefore, a strong 
statement can be made: HMXB/NS microquasars appear as good coun- 
terparts of low-latitude unidentified EGRET sources. 

This statement is followed by some natural questions. 1) May LS 5039 
and LS I +61 303 still contain black holes (BHs)? Although formally 
possible, this is unlikely, because regarding the properties of their X-ray 
emission they should be in the so-called low/hard state (Fender &: Mac- 
carone 2004), but they do not follow the empirical correlation between 
X-ray and radio flux found by Gallo et al. (2003), in the sense that they 
are clearly too radio loud. 2) Why no HMXB microquasars with BHs 
(e.g., Cygnus X-1) are present in the third EGRET catalog? A possi- 
bility is that black hole state changes may play a role, preventing the 
detection of gamma-ray emission when the jet is not present during the 
high/soft state (but this state is very rare in the case of Cygnus X-1 and 
in fact, as discussed by Romero et al. (2002), a low high-energy cutoff of 
a few hundreds in the Lorentz factor of the electrons is compatible with 
the data). 3) Why no Low Mass X-ray Binary (LMXB) microquasars 
are present in the third EGRET catalog? One possibility is that the 
optical companions do not provide the necessary intense UV radiation 
fields needed for an effective IC process to produce high-energy gamma- 
rays. This possibility is strongly model dependent, and not valid if 
Self Synchrotron Compton losses are dominant. The other possibility 
is that LMXBs are in general terms transient objects, that maybe were 
not active during the EGRET viewing periods, while HMXBs tend 
to be persistent systems. 4) Of course, BH and LMXB microquasars 
could still emit high-energy gamma-rays and not be present in the third 
EGRET catalog because of the relatively poor sensitivity threshold of 
the instrument. 

In conclusion, persistent HMXBs containing NSs not experiencing 
state changes are good candidates for the counterparts of the still 
unidentified high-energy gamma-ray sources at low-galactic latitudes 
(approximately up to |6| < 10°, as discussed for LS 5039), and we 
consider that these objects may define a population among UESs. Ob- 
servations with the future missions AGILE and GLAST will confirm or 
reject both the proposed associations between these microquasars and 
the corresponding high-energy gamma-ray sources, and the hypothesis 
discussed above. 
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